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^) • Abstract 

■ We study the competence between chiral symmetry restoration and two flavor color super- 
conductivity (2SC) using a relativistic quark model with covariant nonlocal interactions. We 

■ consider two different nonlocal regulators: a Gaussian regulator and a Lorentzian regulator. We 
find that although the phase diagrams are qualitative similar to those obtained using models 

■ with local interactions, in our case the superconducting gaps at medium values of the chemical 
potential are larger. Consequently, we obtain that in that region the critical temperatures for 
the disappearance of the 2SC phase might be of the order of 100-120 MeV. We also find that for 

I ratios of the quark-quark and quark-antiquark couplings somewhat above the standard value 

. 3/4, the end point and triple point in the T — fi phase diagram meet and a phase where both 

' the chiral and diquark condensates are non-negligible appears. 
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The understanding of the QCD phase diagram has become one of the most interesting issues in the 
' physics of strong interactions. Such phase diagram is relevant to phenomena in the early universe, 
^ in the interior of neutron stars and in relativistic heavy ion collisions. Already in the seventies 

it was suggested that there should be two distinct phases: a low temperature and density phase 
in which quarks and gluons are confined within hadrons and chiral symmetry is broken, and a 
high temperature and density phase (the so-called quark-gluon plasma) in which these particles are 
deconfined and chiral symmetry is restored. Although the possible existence of other phases (e.g. 
the color superconducting phase) was also suggested quite long ago^, this two-phase structure 
became the standard picture of the QCD phase diagram for more than two decades. In recent 
years, however, it has been established that at low temperatures and medium densities several 
other phases might appear (see e.g. Ref. |2j for recent reviews). In particular, for the case of 
two light flavors it has been shown that there should be a non-negligible region in the QCD phase 
diagram where strongly interacting matter is a color superconductor (2SC phase) (2) • Unfortunately, 
due to difficulties in dealing with finite chemical potential, ab initio calculations (as e.g. lattice 
QCD) are not yet able to provide a detailed knowledge of the QCD phase diagram[3]. Thus, most 
theoretical approaches are based on the use of effective models of QCD. Among them, the Nambu- 
Jona-Lasinio modeljSj is one of the most popular. In this model the quark fields interact via local 



four point vertices which are subject to chiral symmetry. If such interaction is strong enough, chiral 
symmetry is spontaneously broken at zero temperature and density, and pseudoscalar Goldstone 
bosons appear. It has been shown by many authors that when the temperature and/or density 
increase, the chiral symmetry is restoredjnj. When effective quark-quark interactions are added to 
the effective lagrangian other phases also appearfTI. As an improvement on the local NJL model, 
some covariant nonlocal extensions have been studied in the last few years0. Nonlocality arises 
naturally in the context of several successful approaches to low-energy quark dynamics as, for 
example, the instanton liquid modeiP and the Schwinger-Dyson resummation techniques [lOJ. It 
has been also argued that nonlocal covariant extensions of the NJL model have several advantages 
over the local scheme. Several studies [TTl 1121 113| have shown that these nonlocal models provide 
a satisfactory description of the hadron properties at zero temperature and density. Recentlvjl4j. 
the characteristics of the chiral phase transition have been investigated within this kind of models. 
The aim of the present work is to extend such studies to the case in which a superconducting phase 
can appear. 

This article is organized as follows. In Sec. II we introduce the model and formalism. Our 
results for some specific nonlocal regulators are presented in Sec. III. In Sec. IV we analyze the 
dependence of our results on the strength of the diquark correlations. Finally, in Sec. V we give 
our conclusions. 



II The formalism 

Let us begin by stating the Euclidean action for the nonlocal chiral quark model in the case of two 
light flavors and SU (2) isospin symmetry. For the purpose of the present study it is enough to only 
consider interactions in the scalar quark-antiquark and quark-quark channels. Thus, one has 



[q{x) qix)]^ - — [q{x) i75r2A2 qc{x)] [qc{x) i-f5r2h q{x)] \ , (1) 

where nic is the current quark mass and T2 and A2 are Pauli and Gell-Mann matrices corresponding 
to the flavor and color groups, respectively. The delocalized quark flelds q{x) are deflned in terms 
of the quark flelds ip{x) as 

q{x) = / d'^y r{x - y) ipiy) (2) 



where the function r(x — y) is a nonlocal regulator. It can be translated into momentum space, 

Lorentz invariance implies that r{p) can only be a function of p^. Hence we will use for the Fourier 
transform of the regulator the form r(p^) from now on. In addition, in Eq.(^ we have used 

qc{x) = 7274 q\x) ; qc{x) = q\x) 7274. (4) 

The partition function for the model is deflned as 
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where the EucHdean action at finite temperature T and chemical potential /i is obtained from 
Eq. by going to momentum space and performing the replacement 



(2^ 



F{p, LOn - in) , 



(6) 



where iVn are the Matsubara frequencies corresponding to fermionic modes, uJn = (2n + l)7rT. As 
in Ref.|l4 we are assuming here that the quark interactions only depend on the temperature and 
chemical potential through the argument of the regulators. To proceed it is convenient to perform 
a standard bosonization of the theory. Thus, we introduce the sigma meson field a and the scalar 
diquark field A and integrate out the quark fields. In what follows we will work within the mean 
field approximation (MFA), in which these bosonic fields are expanded around their vacuum ex- 
pectation values a and A and the corresponding fiuctuations neglected. Within this approximation 
and employing the Nambu-Gorkov formalism, the mean field thermodynamical potential per unit 
volume reads 



T CT^ 1^12 

1 'yMFA " ill 



— I r in 



(2^)3 2' 



(7) 



where Tr stands for the trace over the Dirac, flavor, color and Nambu-Gorkov bispinor indexes. 
The inverse propagator S~^{a,K) is 



S-\a,K) 



(8) 



where 



(9) 



and Vp = r (p ^ -|- — i/x)^). After some straightforward algebra Q'^'"'^ can be more explicitly 



expressed as 



2G ^ 2H 



2T 



d^p 
(2^ 
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{Ap + \Ap\^Y -B, 

rp4 



4C2 



+ In 



^p 



Bp - 4^2 



rp4 



(10) 



where, 

Art — 



^2 + p2^/,2^|Sp|2, 



Bp = 4p'^ + Im^Sp) , Cp = fj, ReSp -|- uon ImSp . (11) 



For flnite values of the current quark mass, ^1'^'"'^ turns out to be divergent. The regularization 
procedure used here amounts to define 



ntMFA 



^MFA _ ^free ^free 
(reg)' 



(12) 



where ^^^^eg) regularized expression for the thermodynamical potential of a free fermion gas, 



free 
(reg) 



-12 T 



d^p 
(2^ 
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In 1 e" 



(13) 
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The mean field values a and A are obtained from the coupled pair of gap equations 

, ^^^ = 0. 



da 

The explicit form of these equations is 




(14) 







^£-J {Ap + \Ap\^f -Bp -A 2{Aj-Bp-iC^) 

oo 

1 - 16 F r 



d^p 
(2^ 



Ap + |App 



{Ap + \Ap\^f-Bp-AC^ 



where, in addition to the definitions given in Ea. pi|) . we have used 



ReSp Re + ImSp Im 



Fp = 2 p'^ ImSp Im Tp 



Gp = fj-Ke Vp + ojn Im 



(15) 



(16) 



It should be noticed that, in general, there might be regions for which there are more than one 
solution for each value of T and /i. In such regions we identify the stable solution by requiring it 
to be an overall minimum of the potential. 

Given the thermodynamic potential the expressions for all other relevant quantities can be 
easily derived. For each flavor the quark-antiquark condensate (V'V') and the quark density pq are 
given by 

m = ^ , p. = -^- (17) 



drur 



dp. 



In the case of the quark-quark condensate an extra source term — ^ '4^c{x)i^b'T2^2'4'{x) has to be 
added to the effective action. It is easy to see that this leads to a thermodynamic potential ^mfa{0 
which has the form given in Ea.H10|) but where Ap has been replaced by Ap — ^. Then, we get 



(18) 



5=0 



Finally, a magnitud which is important to determine the characteristic of the chiral phase transition 
is the chiral susceptibility x- It can be calculated as 



X 



p,2ctMFA 



drrir 



(19) 



III Numerical results for different regulators 

In this section we concentrate on the numerical results obtained for two different regulators often 
used in the literature: the Gaussian regulator and the Lorentzian regulator. In each case we 
have considered G, rric and A as input parameters flxed so as to reproduce the phenomenological 
values of the chiral condensate, pion mass and pion decay constant at vanishing temperature and 
densities jlll I14j . Moreover, we have set H/G = 3/4 as implied by, for example, OGE interactions 
inilZj. The dependence of our results on this ratio will be discussed in the following section. 
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III.l Gaussian regulator 

In this case the regulator is given by 

r(/) = exp(-//2A^)> (20) 

where A plays the role of a cut-off parameter. We have considered two different sets of parameters. 
Set I corresponds to G = 50 GeV-"^, rric = 10.5 MeV and A = 627 MeV while Set II to G = 
30 GeV-^, m = 7.7 MeV and A = 760 MeV. Although for both sets the zero temperature and 
density properties mentioned above are well reproduced, in the case of Set I the quark propagator 
has no purely real poles while for Set II it does. Thus, following Ref.|llj. Set I might be interpreted 
as a confining one, since quarks cannot materialize on-shell in Minkowski space. 

The mean field values a and A as a function of ^ and T are obtained by numerically solving 
Eas. (|15j) . The corresponding results as a function of ^ for various values of T are displayed in Fig^ 
where left panels correspond to Set I while right ones to Set II. It can be seen that, for small values 
of T and // (full lines in Fig^, the system is in the chiral phase for which a 7^ and A = 0. If we 
increase // keeping T fixed, at some critical value of /i there is sudden drop of a and a simultaneous 
sudden increase in A so that we get into the 2SC phase characterized by cj = and A 7^ 0. In 
particular, the values of the diquark gap at the critical chemical potential and T = can be found 
in Table 1. If we repeat the process with a higher value of T something similar happens until we 
reach the "triple point" (3P). For temperatures slightly higher than T3P (dashed lines in Fig^ 
the sudden drop in a and the increase A start to happen at two different values of /i. Between 
these values of /i we have a = and A = 0. Moreover there is no discontinuity in the behavior 
of A as a function of /i. For temperatures above the "end point" (EP) the discontinuity in a 
also disappears (dotted lines in Fig^. Finally for temperatures above the critical temperature for 
/X = 0, Tc(/i = 0), we get a = for all values of /i. In the region corresponding to the crossover the 
transition point is defined by the point at which the chiral susceptibility x has a maximum. The 
positions of the different critical points are summarized in Table 1. Also shown in Fig^ are the 
corresponding quiral (^V) diquark {il^ip) condensates. Their behavior is quite similar to those 
of the mean field values a and A, respectively. It is worthwhile to mention that the fact that the 
chiral condensate approaches some positive value for large values of is due to the substraction 
scheme used to regularize the thermodynamical potential, Ea. ()12() . In fact, it is not difficult to see 
that for finite values of nic the regularized free thermodynamical potential ^j^^^^y Eq. (|13|) . has such 
behavior. 

The corresponding phase diagrams are displayed in Fig|2l Again left panels correspond to Set 
I and right ones to Set II. On the other hand, the upper panels correspond to the phase diagrams 
in the T — fi plane and the lower ones to the diagrams in the T — p/po plane, where the nuclear 
matter density po = 1.3 x 10^ Me^. In all the cases we have indicated with full lines the first order 
transition lines, with dashed lines the second order transition lines and with dotted lines the lines 
corresponding to the crossover between the chiral phase and the weakly interacting quark phase. 

As it is clear from the T — p phase diagrams, at the triple point the three phases can coexist. It 
is interesting to remark that for values of the chemical potential below p^p the transition line (both 
the first order and crossover sections) coincides exactly with that obtained in the absence of diquark 
correlations {H = 0). On the other hand, for p > p^p the first order transition line is different from 
that obtained in Ref.|14j. This is more clearly seen in the corresponding T — p/po diagrams. There 
we have indicated with dash-dotted line the first order transition line corresponding to H = 0. As 
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we observed the existence of diquark correlations increases the size of the mixed phase. Note that 
for T < T3P such mixed phase is composed by the chiral and 2SC phases while for T > T^p it is a 
mixture of the chiral and the free quark gas phases. 

As for the comparison between the results of Set I and II, we see that they are qualitatively 
very similar, with only small quantitative differences in the position of the critical points. 

III. 2 Lorentzian regulator 

The Lorentzian regulator we have considered has the form 



The input parameters are G = 28.38 GeV^"^, nic = 4.57 MeV, A = 940 MeV. Once again, we 
have solved the gap equations for different values of the temperature and chemical potential. The 
qualitative behavior of the mean field values and condensates is very similar to that found using the 
Gaussian regulator (Fig^ and, thus, will not be explicitly shown. Nevertheless it is interesting to 
note that, as shown in Table 1, in this case the T = diquark gap at the critical chemical potential 
Ac is somewhat smaller. The corresponding phase diagrams in the T — fi and T — p/pQ planes are 
displayed in FigEI As in the previous cases, we observe the existence of an "end point" , at which 
the first order chiral phase transition line becomes a crossover line, and a triple point at which the 
three phases co-exist. The position of these points, which are quite similar to those of the Gaussian 
regulator Set II, are given in Table 1. 



IV Dependence of the phase diagrams on the diquark couphng 



In the previous section we have assumed H/G = 3/4 as favored by various effective models of 
quark-quark interactions. However, this value is subject to rather large uncertainties. In fact, 
so far there is no strong phenomenological constrain on the ratio H/G. Thus, it is worthwhile 
to explore the consequences of varying H/G in the range < H/G < 1. Larger values of this 
ratio are quite unlikely to be realized in QCD and might lead to color symmetry breaking in the 
vacuum. Since the results obtained above for different regulators are qualitatively very similar in 
what follows we will only consider the Gaussian regulator with the set of parameters Set II. 

For values of H/G in the range 0.17 < H/G < 0.82 the resulting phase diagrams look qualita- 
tively similar to the one displayed in FigEl although the details (in particular the position of the 
critical points, see below) do depend on H/G. For H/G < 0.17 there is a qualitative change since 
the triple point ceases to exist. In this case, even at very low temperatures, as we increase p, at 
some point we find a first order phase transition between the chiral phase and the free quark gas 
phase. For values slightly above this critical p we have a ~ A = 0. If we continue to increase p 
we find a second order phase transition between the free quark gas phase and the 2SC phase. The 
corresponding T — p phase diagram is shown in the upper panel of Fig0J Note that in this phase 
diagram the transition line between the chiral and the free quark gas phases (both its crossover 
and first order sections) coincides exactly to the one obtained for H = 0. For H/G > 0.82 the 
situation is again qualitatively different, since in such range the triple point and the end point 




(21) 



constant 



6 



merge together. In fact, as H/G comes closer from below to H/G = 0.82 the first order transition 
line that connects the 3P and the EP becomes shorter and shorter and at this particular value it 
disappears. Moreover, for H/G > 0.82 there is small region at low temperatures at which both a 
and A take non-vanishing values. This region is separated from the chiral phase by a second order 
transition line and from the 2SC for a first order transition line. This situation is illustrated in 
the lower panel of Fig0] where the T — fj, phase diagram for H/G = 0.90 is shown. The possible 
existence of this type of phase was already noticed in R,efs. [T5l I16j . 

The behavior of the critical points as function of H/G is displayed in Fig|2I In the upper panel 
we show the position of critical chemical potential /Xc at T = 0. The full line indicates the first 
order while the dashed line the second order one. In the range 0.17 < H/G < 0.82 we have only 
a first order Hc- For values below such fic the system is in the chiral phase while for values above 
is in the 2SC. For values H/G < 0.17 we have that the second order fic, i-e. the point at which the 
second order transition line that separates the 2SC and free quark gas phases meets the /i-axis in 
the T — fi phase diagram, grows rather fast as H/G decreases, signaling the almost disappearance 
of the 2SC phase for very small values of the diquark coupling constant. For H/G > 0.82 a second 
order fic appears again, but now below the first order fic- Thus, between these two critical chemical 
potentials we have a ^ and A 7^ 0. In the lower panel we display the position of the triple and 
end points as functions of H/G. Note that the temperature scale is given to the right while the 
chemical potential scale is given to the left. For values of H/G < 0.17 only the the EP exists. 
Its positions remains independent of the diquark coupling constant up to H/G = 0.82 where it 
meets the 3P, that appears at H/G = 0.19 and which position in temperature (chemical potential) 
increases (decreases) as H/G increases. For values of H/G > 0.82 both critical points transform 
into a single one which position in temperature (chemical potential) increases (decreases) as H/G 
increases. 

V Conclusions 

In this work we have studied the finite temperature and chemical potential behavior of SU{2) f chiral 
quark models with nonlocal covariant separable interactions in both the scalar quark-antiquark and 
quark-quark channels. In our numerical calculations we have considered two types of regulators: 
the Gaussian regulator and the Lorentzian regulator. In all these cases we have set the model 
parameters so as to reproduce the empirical values of the pion mass and decay constant and to get 
a chiral quark condensate in reasonable agreement to that determined from lattice gauge theory 
or QCD sum rules. As for the ratio between quark-quark and quark-antiquark interactions H/G 
which is not well constrained by phenomenology, we have in principle adopted the standard value 
H/G = 3/4 which follows from some models of the QCD interactions. We find that in all cases the 
phase diagram is quite similar. In particular, we obtain that for two light fiavors there are always 
two critical points: a "triple point" at which the second order transition line separating the 2SC and 
normal phases meets the first order transition line which separate the chiral and 2SC phases at low 
temperatures; an "end point" which appears at higher temperatures and at which the first order 
transition line becomes a crossover line. Of course, there is also a critical temperature Tc(// = 0) 
above which the chiral condensate always vanishes. As displayed in Table 1 the values of Tc(^ = 0) 
are in the range 115 — 120 MeV, that is somewhat below the values obtained in modern lattice 
simulations which suggest Tc(^ = 0) ~ 140 — 190 MeVH^ . With this in mind we note that our 
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predictions for the positions of the triple and end points are very similar for all the cases considered. 
Perhaps, the only noticeable difference between the different cases is the prediction for the T = 
diquark gap at the critical chemical potential, where we find values that range from 114 to 182 
MeV. These values are larger than those obtained within models with instantaneous interactions 
j2| . It should be noticed (see Fig^ that in the present case the diquark gap also registers a stronger 
increase with fj, after the phase transition. This leads to rather large values of the gap for chemical 
potentials of the order of 400 MeV, above which strange degrees of freedom have to be taken into 
account. Consequently, for such values of n, we also get larger values (100-120 MeV) for the critical 
temperature needed to go to the free quark gas phase, as it can be easily seen comparing our phase 
diagrams with those obtained in e.g. the NJL model[7l I16j. 

In the final part of this work we have explored the consequences of varying H/G in the range 
< H/G < 1. Given the similarity of the results obtained for the two regulators mentioned above 
we have considered here only the Gaussian regulator with the set of parameters Set II. We found 
that for H/G < 0.17 there is no triple point. On the other hand, for H/G > 0.82 the triple and 
end point merge and a phase where both the chiral and diquark condensates are non-negligible 
appears. It is interesting to remark that the value H/G = 0.82 is quite close to the standard one 
H/G = 0.75 used in most model calculations. Thus, it would be important to sort out possible 
phenomenological consequences of having quark matter with a phase diagram in which the 3P and 
EP coincide. 

In this work we have neglected the strangeness degrees of freedom. To go beyond the values 
of chemical potential considered here and, thus, study for example the Color Flavor Locked phase 
in the context of models with non local interactions their effect have to be included. Work in this 
direction is in progress. 
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Gaussian - Set II 
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Lorentzian 
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114 



Table 1: Critical temperatures, chemical potentials and Ac(T = 0) (all in MeV) for different 
regulators. 
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Figure 1: Behavior of mean fields a, A and condensates {qq), {qq) for the Gaussian regulator, as a 
function of chemical potential for three different values of the temperature. Left panels correspond 
to Set I and right ones to Set II. Full lines correspond to T = 0, dashed lines to T = 67MeV for 
Set I (T = 57MeV for Set II) and dotted lines to T = WOMeV. 
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Figure 2: Phase diagrams for the Gaussian regulator. Left panels correspond to Set I and right 
ones to Set II. Upper panels display the T — ji phase diagrams and the lower ones the T — p/po 
phase diagrams. Full lines indicate first order transition lines, dashed lines correspond to second 
order transition lines and dotted lines to crossover lines. The dash-dotted line in the lower panels 
indicates the section of the transition line corresponding to = 0. 
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to = 0. 
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Figure 4: Phase diagrams for the Gaussian regulator (Set II) and different values of the ratio HjG. 
The upper panel corresponds to the ratio iJ/G = 0.15 while the lower one to B.jG = 0.9. 
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Figure 5: Behavior of the critical points for the Gaussian regulator (Set II) as a function of H/G. 
The upper panel displays the critical chemical potentials at T = 0. The lower panel shows the 
position of the triple and end points. 
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